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SUMMARY 

In this paper the approximation problem for a class of optimal compensators for 
flexible structures is considered. The particular case of a simply supported truss 
with an offset antenna is dealt with. The nonrational positive real optimal 
compensator transfer function is determined, and it is proposed that an 
approximation scheme based on a continued fraction expansion method be 
used. Comparison with the more popular modal expansion technique is 

performed in terms of stability margin and parameters sensitivity of the relative 
approximated closed loop transfer functions. 

INTRODUCTION 

The problem of active stabilization of flexible structures with collocated sen- 
sors/actuators is addressed. In particular, the case of an offset antenna linked by 
a truss to the Shuttle body is considered. A general theory has been already es- 
tablished in [1] as optimal LQG problem for abstract wave equation^. The results 
obtained were applied in [2] to design an optimal compensator for the antenna vi- 
brations suppression after a slewing action, by modelling the truss as a uniform 
Bernouilli beam, simply supported at the Shuttle end, with rate-sensor/actuator 
collocated at the antenna end. The compensator transfer function was determined 
as nonrational, positive real function. This class of' functions was shown to pro- 
vide robust stabilizers for vibrating systems, even in the case of lumped parameter 
systems [3]. Nevertheless, rational approximation schemes are needed in order to 
instrument the compensator. A technique usually adopted is the modal expansion, 
and a typical controller realization can be found in [3] as a bank of filters, centered at 
the frequencies of the system undumped modes. Other methods can be borrowed 
from networks synthesis framework, where the rational approximation of positive 
real functions is a standard problem in telecommunications filters design. Standard 
references on these problems are [5], [6]. Despite a good approximation of some 
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characteristics of the frequency response (amplitude, real part, etc.) can be ob- 
tained, the positive real character of the approximating function is not guaranteed, 
as opposite the modal expansion does for the class of systems considered. This 
is crucial in our control problems since, as mentioned above, positive realness en- 
sures the structure stabilization. Anyway, the main drawback is just the modal 
frequencies computation, obtained by solving a transcendental equation. 

In this paper the compensator transfer function is explicitly computed in the 
vector case of yaw torsion plus roll bending deformation. It is shown to converge 
to a diagonal constant matrix as the control energy increases without bound. Ra- 
tional positive real approximations are deviced via a continued fraction expansion 
technique [7]. Approximations of any order can be easily derived, with coefficients 
straightforwardly related to the system parameters. Moreover, the positive real 
character is guaranteed. The performances of the approximated closed loop trans- 
fer function are evaluated. Comparison with the modal approximation method is 
performed in terms of stability margins and sensitivity to parameters variations. 

THE OPTIMAL COMPENSATOR DESIGN 

We resume in this section the known results about the model and the LQG 
problem for the case of a simply supported uniform Bernouilli beam with an offset 
antenna, with rate-sensor/actuator collocated at the antenna end [2]. Actually, the 
particular case of roll bending deformation (x axis) plus yaw torsion (z axis) is 
considered. In the sequel u^(f,s) denotes the x axis displacement and u v .(t,s) the 
angular displacement about the z axis; t and s 6 (0, /) indicate the time variable 
and the space displacement along the beam axis, respectively. The starting point 
is the following state space model 


Mx(t) + Ax(t)+Du(t) + BN,(t) = 0 (1.1) 


v(t) = B*x(t) + N 0 (t) 


( 1 . 2 ) 


where 

x(f) € Hilbert Space 77 

M : linear bounded, self-adjoint positive definite operator on 77 onto 77, with 
bounded inverse 

A : closed linear operator with domain dense in 77 and range in 77 
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it(-) : the control input G E n , Euclidean n-space 
D : linear mapping E n — ♦ H 

N 3 (-) : white Gaussian noise with spectral density dsl, where I is the identity operator 
of suitable dimension 

N 0 (-) : white Gaussian noise with spectral density d a I y independent of N 9 (.) 
u(.) sensor data. 

The definiton of operators A, M and B, as well as the definition of the space 
H , embody the beam plus antenna dynamical model and the boundary conditions 
stated in [2], For convenience we report these definitions: 



where p is the beam mass density, a the cross sectional area, \ the beam length, 
El#, GI $ are the beam flexural and torsional rigidity respectively, in is the antenna 
mass and r x is the antenna c.o.g displacement, Im is the 2x2 relevant moment of 
inertia matrix of the whole structure. 

We consider the problem of stabilizing the antenna after a slewing action has 
occurred, by determining the control u(.) that minimizes the time average 

i U T ||B-x(i)|| 2 it + A jd IKOir it), A > 0. 

If (A, B) is controllable in [2] it is shown that the optimal compensator transfer 
function is given by 

$0) = apB*(p 2 M + A + jpBB'y'B (1.3) 
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where 

cv — \fd s j d 0 \, 7 VdJd~o T 

'i(f.i) is shown to be a positive real function, thus defining a robust controller [3]. 

THE APPROXIMATION SCHEME 

The compensator transfer function explicit determination. In order to compute 
explicitly the compensator transfer function from eqn. (1.3) consider the following 
expression 


(^M + A + niBB*)x = Bv (2.1) 

which, by taking into account the definitions given in the previous section, can be 
split in the following two relationships 


H 2 M f f + A f f = 0 (2.2) 

f.t 2 Mbb + A b f + ii'yb - v. (2.3) 

By recalling that / = [u*(-) u^(-)]*, (2.2) is solved with the clamped end conditions 
u*(0) = u'^(0) = «,/,( 0) = 0, obtaining 

U^(s) =Ci + C 2 <f>2(s ) 

u^(s) =c 3 sinh (k^s) s €E [0, /), 


where 


<j>i(s) = sinh (hjs) - sin (k^s) f> 2 (s) = cosh(fc^s) - cos (k^s) 

= (paf(EI+)W A) y/iiH k+ = sf^jG n . 

For the constants c* , c 2 the following formulas hold 
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moreover 


«Sr , (-) = [4" , (-) 4 n) (-)i 


>i(0 0(0' 

-1 


-0l(O 02(0- 


> 2 ( 0 . 


n = 0 , 1,2 


Now it can be shown that A^f can be actually expressed as T(f.i)b [4] so that 
( 2 . 1 ) becomes 


{ixMb + T(f.t) + ^I)b = v 

and the compensator transfer function is defined as 

'J'(m) = otn(fi 2 M b + T(n) + (2.4) 

where I is the identity operator. Recalling the definition of a and 7 , we see that 

lim 4/(p) = d 3 / d 0 I 

so that we have the “direct connection” if the control energy increases without 
bound. 

We list below the non zero entries of the 3x3 matrix T(/j) 


rr f..\_ VT ,.3 s mh(^0 cos (^0 + sin(^>0 cosh (40 


TnW =TM = El.kl 


smh(k ( f,l)cos(k^l) — sin(A:^/)cosh(fc^f) 
i-2 2U«) -EL^k^ 


T 33 (i.i) =GI^k t /,coth(k^l) 

A = — 1 + cos\i(k^l)cos(k,f,l). 
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T(p) is a meroinorpliic function; i.e. it is analytic in all the plane but in a countable 
set of points, where it has polar singularities [8]. 

The continued fraction expansion approximation. In order to approximate (2.4) 
note that if a matrix is positive real so is its inverse; consequently it is of great 
simplification to work on 


-(p. 2 M b + T(^) + 7 piI). (3.1) 

Next, since for p = the function T(p) is real, we see that the real part of (2.5) 
does not depend on T(/u), so that whatever approximation we device, a positive real 
approximate function is obtained. Thus let us concentrate on T(fx). 

The approximation scheme proposed consists of a continued fraction expansion 
of a meromorphic function f(z ), i.e. 


f(z) = r 0 (z) + 


r 2 ( z ) + 


r 3 (z) + 


r 4 (z) + . . . 


(3.2) 


where the 7 ^( 2 ), i = 0, 1,2, ... are rational functions of finite degree suitably defined, 
and are called “convergents” of the continued fraction. The meaning of (3.1) is the 
following [7]: denoting f n (z ) the function obtained by considering n convergents on 
the r.h.s we have 


lim f„(z ) = f(z) 

n— mx) 

for every value of 2 in the complex plane. 

The functions 7 ^( 2 ) can be determined according to the following algorithm. 

step 1. Consider the Laurent expansion about the origin of f{z ) (see e.g. any standard 
book on complex analysis or passive networks synthesis) 

no 00 

f(z) = Yl a -kZ~ k +5^0*2*. 

Jfc=0 k = 1 
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The finite sum is called “singular part” of f(z ) at the origin, n 0 is the multi- 
plicity of the pole at the origin and set r 0 (z) = a -kZ~ k ] it can eventually 

be a constant (n 0 = 0); i.e. the function is regular at z = 0. 
step 2. Compute 

/i 0 ) = /O) - '0(2), 

then fi(z ) has the same poles as f{z ) but the pole at z = 0. Moreover it holds 

f(z) =r 0 (z) + fi(-) (*) 

lim fi(z) = 0 . (it) 

z — ►O 

Note that (i) is not a local expansion but an exact representation of f(z) with 
its singular behaviour at z — 0 explicited. 
step 3, Consider 

/2O) = 1//1O); 

it has a pole in the origin according to (it), then use the argument of steps 1 
and 2 to determine / 2 (z) = ^i(z) + fa(z) with r\{z) = Sjt=o ^-k z ~ k and ni is 
the multiplicity of the pole at z = 0 of / 2 (z ). Since /i(^) = 1/ / 2 (z) we have 

2) = r 0 (*) + 7 . * . ( z . 

ri{z) + / 3 (z) 

It is clear at this point how to obtain expansion (3.1) by repeating step 3 and 
determining recursively all the convergents r*(z), k = 0,1,2.... For completeness 
let us see how to easily obtain the coefficients of the singular part of the Laurent 
expansion at the origin of a meromorphic function f(z) 

a_ n = lim z n f(z) 

2 — 0 

= hm z k ( f(z) 
z — 0 l ' 

The algorithm described can be applied to each entry of T(j.i) obtaining an 
expansion of the following type 


n-k 


n 1+1 )• 


t — 1 


R iJ ( a0 


K' J + 


* 1 ? 

h 4 


+ k'' } 


10 


(3.3) 


+ 


k' J 

221 1 

a ^ 20 


+ 


k' }} 


+ k 


30 


1 

+ — 


489 



and consequently the n-th order approximation T n (f.i ) of T(/i) is obtained as 


TM 


. o o 


Here R'^(ii) is obtained by (3.3) including n convergents in the continued fraction. 
This number needs not be the same for all the entries, and in this case the order 
for T„(/i) is determined by the highest value of n used. 

In the following table are reported the coefficients of the continued fraction 
expansion up to the 4-tlx order of the entries of T(p), for a particular choice of the 
system parameters (see [9]). 

Table 1. Coefficients of the Continued Fraction Expan- 
sion of entries. 



R 1 ’ 1 

R 1 ’ 2 

R 2,2 


ko 

12 

-6 

4 

1 

k n 

-2.7 

19.1 

-105 

3 

&10 

.0020 

-.028 

.18 

.2 

^21 

411815.38 

-74697.1 

-10.2 

-175 

^20 

-451.03 

50.57 

21610.7 

-7.7 

^31 

94.35 

714.4 

-3434.97 

35.04 

&30 

.0185 

-.104 

.4 

.01 

k 4i 

5501900 

-921286.70 

244973.48 

-782.41 

/c 40 

321.53 

44.58 

-10.15 

-7.08 


We stress that the coefficients of the approximations are easily obtained just 
by computing limits in the origin of suitably defined functions, and are simple com- 
binations of the system parameters which appear in the coefficients of the function 
to be approximated. As a result we have a procedure which is numerically robust 
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since no solution of transcendental equations are required, no derivatives determina- 
tion, no inner products computations, no variable transformations as usually occurs 
in the most popular approximation methods. This straightforward relation to the 
system parameters determines also a low sensitivity to parameters variation of the 
approximated closed loop transfer functions, as we will show next. 

For convenience, we report briefly the modal expansion technique (see [3]). 
Let u> k, k = 1,2, ... , the system undamped modes, and ij’k, k = 1,2,..., the M- 
orthogonal eigenvectors, normalized as [Mipki^k] = 1. The following compensator 
transfer function modal approximation can be deviced 



v 2 +^l + mhk' 


Re/i > 0, 


where 

b k = B*xl> k , b k k = 

Thus we have a bank of band pass filters centered at the undamped inodes. 

CLOSED LOOP TRANSFER FUNCTION PERFORMANCES EVALUATION 

In this section we compare the performances of the approximation method pro- 
posed with respect to the modal expansion technique, widely used in this field. 
In particular the stability margins of the approximate closed loop transfer func- 
tion are considered and the sensitivity of this performance index toward system 
parameter variation is evaluated. 

As it is well known the stability properties mentioned can be derived by exam- 
ining the frequency behaviour of the following function 

S(w) = det(J + P(w)*(u;)), 

where P(w) is the system transfer function, in our case defined as 

P(u;) = (T(u>) — u> 2 Mb)~ l . 

Actually, we are more interested in the sensitivity of stability margins with respect 
to parameters variations. In Fig. 1 amplitude and phase plots of the diagonal 
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entries of F(u>) = P(u;)'I>( u>) are reported for a second order continued fraction 
approximation of 'f'(o’). In Fig. 2 are reported the corresponding plots for a modal 
approximation based on two band pass filters centered on the first two modes. 
Particularly interesting are the phase plots showing that in the first case we obtain 
a higher phase margin, practically equal to n/2 for all frequencies. A similar result 
holds for the other entries of F(u>), thus giving a description of the stability features 
achieved in both approximation schemes, avoiding to get through the complexity of 
the function S(uj). 

Moreover, in Fig. 3 are reported the plots of the sensitivity of the phase functions 
considered, with respect to the variation of the parameter 9 = -gj^. Here we note 
that the continued fraction approximation shows a better performance in terms of 
robustness than the modal expansion. 

CONCLUSIONS 

The continued fraction method proposed allows to approximate any meromor- 
phic function by operating simple computations on the coefficients, i.e. the deter- 
minations of the limit in the origin of suitable functions derived by the assigned 
one. This results in a good performance of the approximation in terms of stability 
margins and robustness of the approximate closed loop transfer function, lliis fea- 
ture is highlighted by comparing the mentioned characteristics with the analogous 
one obtained by using the more popular modal approximation scheme. 
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